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ABSTRACT 

The evolution of the quasar luminosity function is one of the basic cosmological 
measures providing insight into structure formation and mass assembly in the Universe. 
We have conducted a spectroscopic survey to find faint quasars (—26.0 < M1450 < 
—22.0) at redshifts z = 3.8 — 5.2 in order to measure the faint end of the quasar 
luminosity function at these early times. Using available optical imaging data from 
portions of the NOAO Deep Wide-Field Survey and the Deep Lens Survey, we have 
color-selected quasar candidates in a total area of 3.76 deg^. Thirty candidates have 
R < 23 mags. We conducted spectroscopic followup for 28 of our candidates and found 
23 QSOs, 21 of which are reported here for the first time, in the 3.74 < z < 5.06 redshift 
range. We estimate our survey completeness through detailed Monte Carlo simulations 
and derive the first measurement of the density of quasars in this magnitude and redshift 
interval. We find that the binned luminosity function is somewhat affected by the K- 
correction used to compute the rest-frame absolute magnitude at 1450A. Considering 
only our i? < 23 sample, the best-fit single power-law (<I> oc L^) gives a faint-end 
slope P = —1.6 lb 0.2. If we consider our larger, but highly incomplete sample going one 
magnitude fainter, we measure a steeper faint-end slope —2 < (3 < —2.5. In all cases, we 
consistently find faint-end slopes that are steeper than expected based on measurements 
at z ~ 3. We combine our sample with bright quasars from the Sloan Digital Sky Survey 
to derive parameters for a double-power-law luminosity function. Our best fit finds a 
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bright-end slope, a = —2.4 it 0.2, and faint-end slope, /? = —2.3 it 0.2, without a well- 
constrained break luminosity. This is effectively a single power-law, with f3 = — 2.7ib0.1. 
We use these results to place limits on the amount of ultraviolet radiation produced by 
quasars and find that quasars are able to ionize the intergalactic medium at these 
redshifts. 



Introduction 



Understanding the evolution of quasars has been a subject of great importance since their 
discovery over four decades ago. In particular, there is now substantial, multifaceted, and growing 
evidence for a correlati on between the forma t ion and evolution of gala, xies and their central su- 
permassi ve black h oles (jMagorrian et al.lll998l: Ferrarese &: Merrittll2000l : also, see the proceedings 
edited by|Ho||200J, or the review bv lDiorgovskill2005l . and references therein). Studies of the evolv- 
ing QSO populations are important not only for their own sake, but also for providing insights into 
the formation and evolution of massive galaxies in general. 

While quasars represent a relatively minor contributor (~ 10%) to the overall energetics of 
the post-recombination universe, they dominate at high energies, and their radiative and mechan- 
i cal feedback may significantly affect the formation and growth of their hosts and companions 
(jSilk &: Reed Il998l : IWvithe fc Loebl l2003bl ) . Quasars are a significant contributor to the meta- 
galactic ionizing radiation field at any redshift, although their role relative to that of star- forming 
galaxies has likely changed over the history of the Universe. Quasars might still be important 
producers of the metagalctic UV radiation at high redshifts (z > 4 — 5) as one approaches the end 
of the reionization era. 

Although we observe a decline in the density of quasars by a factor of ^ 40 between the peak of 
quasar activity at z ~ 2.5 and the end of reionization at z ~ 6 (jFan et al.ll2001al ). this evolution has 
been traced only by the most luminous sources at redshifts beyond z ~ 3. The primary observable 
which constrains the evolution of quasar populations and their effects on their environment is the 
quasar luminosity function (QLF) as a function of redshift. At low reds hifts, the QLF is well- 
repre sented by a double power- law, $(L) = $./L.[fL/LJ" + fL/LJ^]-^ (lBovle et al.lll988l : [pdl 
19951 ). The bright-end slope, a, has typical measu red values around —3.4 at redshifts z < 2.5, and 



flattens toa~— 2.5at2;~4 (|Fan et al.ll2001bl ). T he faint-end slope , f3, is typically measured 
to be around —1.7 at low redshifts (z < 1), although iHunt et al.l (|2004l ) find a shallower value of 
P « -1.2 at z ~ 3. 

There has been considerable progress recently in the theoretical understanding of the shape 
of the QLF and its evolution. The emerging picture is a complex interplay of QSO lifetimes and 
luminosity histories, powered by dissipative mergers, with the bright end of the QLF dominated by 
short-lived phases of QSOs radiating near the peak of their l uminosities and the faint end determined 
by the distribution of QSO fueling lifetimes and feedback ( Hopkins et al. 2005 . 2006al jbl). Thus, in 
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this model the faint end of the QLF at high redshifts represents a critical observational constraint 
on the early formation history of massive black holes, their contribution to the reionization, and 
feedback processes affecting the formation of their host galaxies. 

The true shape of the QLF at z > 4 is currently unknown, largely due to the flux limits of most 
large-area surveys to date. The recent availability of deep, relatively wide-field, multicolor, optical 
surveys in the public domain have enabled the search for faint quasars at these high redshifts and 
the de termination of the fa, int end of their QLF. In this paper, we utilize the De ep Lens Survey 



(DLS: IWittman et al.l 12002 ) and the NOAO Deep Wide-Field Survey (NDWFS; Ijannuzi k Dev 



19991 ) which go ~ 4 magnitudes deeper than the Sloan Digital Sky Survey (SDSS) and other large 



area sky surveys, to measure the faint end of the QLF at z ~ 4. 

The organization of the paper is as follows: in Section 2, we review the NDWFS and DLS 
survey data that were used for photometric candidate selection. Section 3 presents the simulations 
of quasar colors we have performed to obtain the QSO candidate color selection criteria. In Section 
4, we describe the followup spectroscopy observations. In Section 5 we compute the QLF and 
analyze its shape in conjunction with published measurements at brighter magnitudes. In Section 
6 we examine the implications of our QLF by computing the contribution of quasars to the UV 
radiation field at 2; ~ 4. We discuss our results in Section 7. We use standard cosmological 
parameters throughout the paper: Hq = 70 km s~^ Mpc~^, JIm = 0.30, and Qa = 0.70. 



Candidate Selection 



We select our quasar candidates using the standard technique of looking for objects whose 
colors are outliers from the stellar locus. At 2: ~ 4, o ur redshift of interest, -B , R, and / (or 



effectively and efficiently separate quasars from stars (jKennefick et al 



1995 



Richards et al 



2OO2I ). We make use of publicly available deep imaging data from the NDWFS and the DLS, which 
we describe below. 



2.1. NOAO Deep Wide-Field Survey 

NDWFS is a deep imaging survey of two 9.3 deg^ fields, Bootes and Cetus, in three optical 
and infrared bands {Byy, R, I and J, H, K). For our candidate selection, we obtained imaging 
data in the Bootes field from the Third Data Release (DR3) which is publicly available from the 
NOAO Science ArchivtEl. 

The custom B^r filter has high transmission over most of the wavelength range commonly 
covered by U and B filters, as can be seen in Figure [H Optical imaging for the data available 



|http : //www . eir chive . noao . edu/ndwf s/ 



- 4 - 



in NDWFS DR3 was obtained with the 4 meter Kitt Peak National Observatory (KPNO) Mayall 
telescope and MOSAIC-I Wide Field Imager. The Bootes field is split into 27 partially overlapping 
subfields, each roughly 36'x36' on the sky, with imaging products delivered on a common (tangent- 
projected) scale of 0'.'258 pixel"^. The survey design aimed for 5a point source detection limits of 
Bw ~ 26.6, R « 26.0 and / « 26.0 (AB mags). 

We use optical {Bw, R and /) imaging for five of the 27 NDWFS subfields, selected to optimize 
the seeing conditions and total exposure time in all three filters. The five subfields cover a total 
area of 1.71 deg^ and their details are summarized in Table [TJ 

The B\Y^ R, and / images for a given subfield have different sizes. We registered the images 
and trimmed them to their common overlapping area (lis ted in column 4 of Tab le [1]) . We then 



extracted source catalogs using SExtractor version 2.4.4 (iBertin Sz ArnoutsI 119961 ) in dual image 
mode with i?-band as the detection image and forcing the measurements in Bw- and /-bands. For 
detection in SExtractor, the i?-band image was smoothed with a 3 pixel wide Gaussian filter, and 
the detection threshold was set to 2a of the sky noise. 

Photometry was performed in 3" diameter apertures. To compensate for differences in seeing 
between the three bands (see Tabled]), this aperture magnitude was corrected to "total" magnitude 
for each band separately by subtracting an aperture correction. To determine the aperture correc- 
tion in each image, we compared the isophotal magnitude for >200 bright stellar objects to their 
aperture magnitude and measured the median difference. This correction is typically < 0.1 mag 
in e very band. We also c orrect all the magnitudes for interstellar extinction using the dust map 



from lSchlegel et al.l (|l998l ). For our selection, we work with these extinction and aperture-corrected 
magnitudes in the ABi, system, converting the published Vega-based NDWFS zeropoints to ABiy 
magnitudes by adding 0.0, 0.26 and 0.53 mag to the B]y, R and / zeropoints, respectively. 



2.2. Deep Lens Survey 

DLS is closely related to the NDWFS in that the optical imaging has been obtained on the 
same 4 meter Mayall Telescope (at KPNO) and Blanco Telescope (at CTIO) using the same wide- 
field optical imagers Mosaic-I and Mosaic-II as the NDWFS. DLS imaging will ultimately cover a 
total of 28 deg^ in four optical bands {B, V, R and z'), separated on the sky into seven fields of 
4 deg^ (labeled Fl - F7). Each of these fields in turn is divided into 9 subfields, or "pointings", of 
35.1' X 35.1'. The final image products released by the DLS are already registered in pixel space. 
Since the scale for these data is also 0'.'258 pixel"^, the final images have sizes 8192 x 8192 pixels. 
The DLS observational strategy was to obtain the i?-band images on the best of the observing 
nights, while the B, V and z' data were taken in conditions of mixed relative quality. In addition, 
the survey was designed to reach 12,000 seconds (20 x 600 seconds) in the B, V, and z' filters and 
18,000 seconds (20 x 900 seconds) in the iZ-band. 

We used the images from the third public data release of DLS data for 6 pointings, covering a 



- 5 - 



total area of 2.05 deg^. The imaging data of the NDWFS and DLS differ in principle only in the 
filters used and fields surveyed; therefore, we applied the exact same procedure and parameters for 
catalog extraction as for the NDWFS data, including thresholds, apertures, and subsequent seeing 
and extinction corrections (see §2.1). The positions and several properties of the field pointings 
used in our survey are presented in Table [2l 

We extracted catalogs using SExtractor in dual image mode, with the i?-band images as the 
detection image. Photometric measurements were forced on the same pixel positions in the B, 
V and z' images. The DLS data provide photometric zero-points. The B, V, and R images are 
calibrated to the Vega system, while the z'-band images are calibrated to the SDSS AB^^ system. 
We convert from Vega magnitudes to ABi, by adding —0.09, 0.0, and 0.2 to the B, V, and R zero- 
points, re spectively. We also c orrect all the magnitudes for interstellar extinction using the dust 
map from lSchlegel et al.l ( 1998 ). Although our main selection filters are B, R and z' , we generated 
catalogs for the F-band data to further discriminate between our candidates and increase our 
efficiency, as we describe below. 



2.3. Field-to-field Consistency 

We examined the position of the stellar locus in the {R—z') vs. (B—R) and {R—I) vs. {B^—R) 
color-color space for the DLS and NDWFS data, respectively, to verify the consistency of the zero- 
points. We found discrepancies between the fields that, while small, are larger than the expected 
photometric errors. 

To determine the zero-point shift required to bring our photometry onto a standard system w e 



invoke a method used in the photometric quality assessment of SDSS. Following llvezic et al.l (|2004l ). 
we defined a principle axis in the aforementioned color-diagrams that follows the main sequence by 
iteratively fitting a line through the observed locus of well-detected stellar objects in each field. We 
determined the origin of the stellar locus from the drop in density of stars along the determined 
track. Figure [2] shows the the positions of the main sequence color tracks for the five NDWFS fields 
on the left and the six DLS fields in our survey on the right. The dashed lines indicate the principal 
axes, tracking that stellar locus, and the circles represent the origins of the locus. 

Two DLS fields, F2p31 and F2p33, have bluer {R — z') colors, and F2p31 is redder than the 
other fields in (B — R). In addition, the photometric quality of F2p31 may be slightly poorer 
that the other five DLS fields as the stellar locus has an rms of ~ 6% compared with 3% for the 
remaining fields. We apply a correction to the z'-band zero-points of —0.15 and —0.18 to F2p31 
and F2p33, respectively. We also apply a z'-band zero-point shift of —0.05 magnitudes to F2p31. 
The NDWFS data were consistent within photometric errors; therefore no correction was applied 
to those fields. 
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3. Color Selection Criteria Using Simulated Quasar Spectra 

To define color selection criteria for quasar candidates at z > 4 and determi ne our sele ction 
function, we sirnulate d quasar spectra following the general procedure outlined in iFanI (jl999l ) and 
Richards et alj (j2006l ) and calculated their colors for the NDWFS and DLS filter combinations. 
We generated a library of simulated quasar spectra, starting with a quasar template spectrum 
constructed fro m a combina ti on of the Hubble Space Telescope (HST) radio-quiet quasar composite 
spectrum from iTelfer et al.1 (|2002l ) which covers 300 - 246lA in the quasar rest frame and the 
SDSS quasar composite spectrum (jVanden Berk et alJ 120011 ) spanning rest-frame 800 - 9200A. 
This combined template uses the HST composite below 1275A and the SDSS composite above 
2000A. The spectra are averaged in the overlap region of 1275 - 2000^0. 

We simulated the attenuation of each model spectrum by intervening neutral hydrogen ab- 
sorbers by shifting it to a given redshift and cr eating random r e alizat ions of a population of discrete 
absorbers following the procedure outlined in iBershady et alJ (|l999l ). The absorbers attenuate the 
template quasar spectrum with their Lyman series absorption lines and the Lyman limit break. 
Two separate distributions of absorbers are used, divided up according to their column density 
(Nhi), corresponding to "weak" absorbers and Lyman Limit Systems. We did not attempt to sim- 
ulate the effect of any dust reddening or absorption by neutral hydrogen intrinsic to the quasar 
host galaxy. 

We account for the intrinsic diversity of quasar spectral properties by varying the power- 
law slope of the continuum and the strength of the Lya emission line. At UV-optical rest-frame 
wavelengths, the QSO continuum is best described by a broken power-law. Shortward of Lya, we 
assign the spectrum a power-law index drawn fr om a Gaussian dis tribution centered on ai = —1.57 
with a standard deviation, cji, of 0.17 {f,y oc v"'; iTelfer et al.ll2002l ). Longward of Lya we assign the 
continuum a po wer-law index drawn from another Gaussian distribution centered on a2 = —0.5 
with (72 = 0.30 (jVanden Berk et al.ll200ll : iRichards et al.ll2003l ). We assume that the values of the 



two spectral slopes are uncorrelated. In addition, we vary the eq uivalent widt h, EW, of the Lya 
emission line assuming a Gaussian distribution of EW=65 it 34A ( Wilkeslll986l ). and again assume 
that it is uncorrelated with any of the other spectral properties. 

At each Az = 0.1 in the redshift range 3.5 < z < 5.2, we created 200 simulated quasar 
spectra, each with a single drawing from distributions of the two power-law continuum slopes, the 
Lya emission line EW, and one realization of neutral hydrogen intervening absorption systems. 
Convolving these spectra with the filter responses for the NDWFS and DLS filters {Bw, R, I and 
B, V, R, z', where the two R filters are identical) we obtain simulated quasar "color tracks" shown 
by the dots in Figure O 



Based on the colors of our model quasar spectra we adopt a set of color cuts for each of our 



^This combination is available at 



littp://www.pha.jliu.edu/~rtl9/composite/"] 
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datasets, NDWFS and DLS. These are shown with dashed hnes in Figure [3j For the NDWFS 
survey, we impose the foUowing selection criteria: 

{Bw -R)> 2.2 (1) 

and 

(R-I) < 0.25 X {Bw -R)- 0.3. (2) 
For the DLS survey, we use the following color cuts: 

{B-R)> 2.2 (3) 

and 

{R- z) <0A5 X {B - R) -0.59. (4) 

The DLS survey provides a fourth filter, V, which we use to impose an additional color criterion, 
based on the colors of the aforementioned simulations: 

{B-V)> 1.0. (5) 



3.1. Morphological Criteria 

In addition to having the right colors, we required our quasar candidates to be point sources. 
We use the full width at half-maximum (FWHM) to discriminate between stellar and extended 
sources. The left panel of Figure H] shows the FWHM as measured from the i?-band image in 
SExtractor as a function of i?-magnitude. A locus of stellar sources with FWHM ~ 1" is obvious. 

We sought to determine a suitable morphological cut of FWHM that would eliminate galaxy 
contamination without sacrificing completeness. To do this, we used the mkobjects task in the 
IRAF artdata package to insert 4200 objects with stellar profiles and 17.75 < R < 24.50 into each 
field image. We re-extract our catalogs using the same SExtractor parameters, and examine the 
FWHM distribution of our simulated point sources. This distribution is shown in the right-hand 
panel of Figure H] for the same field. The solid-line histogram shows the distribution of FWHM for 
simulated point sources, while the dotted line shows the FWHM distribution for sources brighter 
than R = 22 and the dashed line shows the FWHM distribution for sources with 22 < R < 24. 
For each field, we choose as our cutoff the FWHM that retains 80% of our simulated sources. This 
cutoff, calculated to be 1.145" for NDWFS J1426+3236, is shown as horizontal dashed line in the 
left panel of Figure S] and as a vertical solid line in the right-hand panel. 



3.2. The Final Candidate List 

At z > 4, neutral hydrogen absorption from the Lya forest often causes B-hand drop outs, 
even at the deep sensitivities of NDWFS and DLS. Creating our catalogs in SExtractor's dual- 
image mode ensures that we measure a flux for each object that has an i?-band detection, even if 
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an object's flux drops below the B\y, B or V detection limits. Nevertheless, the magnitude assigned 
to a drop-out object may be unreliable. 

We define the limiting magnitude of each field based on the standard deviation per pixel for 
a representative region in each image. We compute the uncertainty in total flux in a 3" aperture, 
with area A, where each pixel has lapix counts: fia = apixVA. We compute the limiting magnitude 
from this flux for the Bw fields in NDWFS and the B and V field in DLS. We assign these limits 
to any object whose SExtractor-determined magnitude is fainter than these limits. We then use 
these limits when applying the color cuts in our candidate selection. We note that by ignoring the 
fact that pixel values are correlated in these images we may be underestimating the true noise. 
Furthermore, this estimate does not take into account large scale variations from the flat field and 
other systematic effects. Nevertheless, this value serves as a practical lower limit on the magnitude 
in a band without a detection. Since the limiting magnitudes that we compute are B, B\y > 27 
mags and our color criteria require that B, Bw — R > 2.2, our survey limit of i? = 24 ensures that 
all undetected objects meeting our color criteria will be selected. 

We ran our catalogs through these selection criteria and ended up with 80 quasar candidates 
in the five NDWFS fields and 74 candidates in the six DLS fields, strictly based on their cataloged 
properties. We then examined image cutouts of the objects and eliminated three candidates from 
the NDWFS sample and three candidates from the DLS sample that were image artifacts mistaken 
as objects by SExtractor. We are therefore left with 77 (16 with R < 23) NDWFS candidates and 
71 DLS (14 with R < 23) candidates. We plot these objects as circles in Figure [3], indicating with 
a right-pointing arrow B\y and S-band drop outs. We list the candidates with i2 < 23 in Tables [3] 
and [Hand make the full candidate list available in electronic table format. 



4. Spectroscopy Results 

We obtained 28 spectra for the DLS and NDWFS candidates and found 23 high-redshift quasars 
(the remaining objects were three Lyman-break galaxies at z ~ 3.9, one Type II AGN at z ^ 0.2 
and one featureless, unidentified spectru m). All but two of our spectra were taken with the Low- 



Resolution Imaging Spectrometer (LRIS: lQke et al.lll995l ) on the Keck I telescope on Mauna Kea 



Five quasars were discovered on UT 2005 November 28 and 29. Only the red camera was used, with 
the 400 lines mm~^ grating blazed at 8500 A. Sixteen additional quasars were discovered during a 
run on UT 2006 May 20 to 22. Again, only the red camera was used, this time with the 600 lines 
mm^^ grating blazed at 7500 A. The spectra were obtained one at a time using the long-slit mode 
at the parallactic angle. 

We reduced the LRIS spectra using BOGUS, an IRAF software package developed by Stern, 
Bunker, &: Stanforcifl to reduce LRIS slitmask data. We made minor modifications to the code 



^Available at |http : / /zwolf kinder . jpl . nasa . gov/~ stern/homepage/bogus . html 
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to accommodate the single slit mode of our observations. BOGUS reduces the two-dimensional 
spectra, applying a bias-correction, flat-fielding, cosmic ray removal, and sky subtraction. Most 
of our spectra were observed in one exposure at a single slit position. A handful of spectra were 
dithered in two positions along the slit. For these objects, BOGUS applied a fringe correction as 
well. The shain_r routine (part of the BOGUS package) extracts the one dimensional spectrum and 
performs wavelength and flux calibrations. Since the fringing can be quite strong at the position 
of sky lines at red wavelengths, we interpolated over strong sky lines which appeared in our final 



spectra using the skyinterp IRAF tasl 

Two additional spectra were taken on UT 2006 March 25 and 26 with the Multi-Aperture 



Red Spectrometer (MARS; iBarden et al.l l200d ) on the Mayall 4 meter telescope at Kitt Peak. 
MARS is an optical spectrograph that uses a high resistivity, p-channel Lawrence Berkeley National 
Laboratory CCD with little fringing and very high throughput at long wavelengths. We obtained 
spectra of DLS J105227.9-054234 and DLS J105348.8-053319 using the 1'.'7 wide long slit, OG550 
order-sorting filter, and the VG8050 grism. Across much of the optical window this instrument 
configuration delivers resolution R ~ 1100 spectra, as measured from sky lines filling the slit. The 
data were processed following standard optical slit spectroscopy procedures. The nights were not 
photometric, b ut relative flux calibration was achieved with observations of spectrophotometric 



standards from iMassev &: Gronwalll (jl990l ) obtained during the same observing run. 



Figures [6] and [7| plot the 23 quasars discovered by our survey, sorted by decreasing redshift. 
We mark with a vertical dotted-line expected positions of prominent ultraviolet quasar emission 
lines: Lya 1216, N V 1240, Si IV 1400, C IV 1550, and C III] 1909. 



4.1. NDWFS 1433-^3408 and DLS 1053-0528 



These quasars, at z = 3.88 and z = 4.02, respectively, show prominent N IV] A1486, an 
extremely rare emission line. When seen in QSOs, this line is normally accompanied by strong N V 
A 1240 and is a consequence of high nitroger i abundances. Such objects are dubbed "ii i trogen-rich" 
and make up at most .2% - 0.7% of QSOs dOsmer fc Smithlll98d : [Baldwin et aDl2003l : Isentz et al. 



20041 : 1 Jiang et al.ll2008l ). The QSOs in our survey, however, do not seem to fit into this population. 
Their N IV] A1486 equivalent widths (EWs) are orders of magnitude higher than nitrogen-rich 
QS Os, and their correspo nding N V A1240 is weak, and may even be absorbed in DLS 1053—0528 
(see lGlikman et al.ll2007l . for more details on these objects). 



Strong N IV] A1486 emission has been seen in the spectrum of the Lynx arc, a gravitationally 
lensed H II galaxy at z = 3.357 (iFosbury et al.l 120031 ). Here this line can be explained by a hot 



''SKYINTERP is part of the WMKONSPEC package for reducing long-slit near-infrared 
spectra, where imperfect sky subtraction can often produce spectra with strong sky lines 
( |http: //www2 .keck.hawaii . edu/inst/nirspec/wmkonspec/index .htmlj ). 
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(80,000 K) blackbody caused by a top-heavy initial mass function (IMF) as the source of ionizing 
flux - suggestive of early, metal-poor star formation, e.g., Population III stars. 



5. Computing the Luminosity Function 

The luminosity function (LF) is defined as the co- moving volume density of objects as a 
function of luminosity. In the case of the quasar luminosity function (QLF) it is customary to use 
the absolute magnitude at 1450A, M1450. To construct the QLF for the quasars in our survey, we 
must first account for the incompleteness of our selection and spectroscopic followup. We must also 
compute -/Vfi45o for each of our quasars. We describe these steps below and compute the QLF for 
our survey. Table [5] summarizes our results from the following sections. 



5.1. Survey Completeness 

To compute the quasar luminosity function, we must be able to correct our quasar counts for 
incompleteness in our selection technique. To determine the completeness of our survey as a function 
of redshift and R magnitude, we follow what is now a stan dard methodology of simulating the 



selection probability of quas ars in the NDWFS and DLS fields (jWarren et al.lll994l : iKennefick et al 



1995bl : [Richards et al.ll200(Tl ^ 



We begin with our simulated quasar spectra, described in Section [3l in a redshift vs. R- 
magnitude grid. The redshift grid points are separated by 0.1 in the range 3.5 < z < 5.2, and 
the magnitude grid points are separated by 0.25 magnitudes in the range 17.75 < R < 24.50. At 
each grid-point we simulate 200 spectra generated from a random drawing of distributions of the 
two power-law slopes, Lya emission line EW and one realization of neutral hydrogen intervening 
absorption systems. The final set of simulated quasars contained 100,800 spectra that evenly 
samples our parameter space. These spectra are convolved with the NDWFS and DLS filter curves 
to produce tables of magnitudes for all the filters in a given survey, as functions of i2-magnitude 
and redshift. 

The NDWFS and DLS images have ~ 45, 000 and ~ 60, 000 - 70, 000 sources detected, re- 
spectively, in each of their fields. We randomly distributed the 100,800 simulated quasars into 
twenty-four subsets of 4200 quasars, to keep the number of added sources under 10% and avoid 
over crowding. To properly insert the simulated quasars into our images as point sources, we deter- 
mined the shape of the point spread function (PSF) by measuring several hundred well-detected, 
unsaturated, stars in each image using the IRAF task imexamine and determining the cr-clipped 
median values of a and (3 in the Moffat profile model, 

/ = 4(1 + (r/af)-^, (6) 



where Ic is the peak value of the intensity, /. We then inserted 4200 quasars into the NDWFS 
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and DLS images, using the IRAF mkobjects task in the artdata package, in all filters at the 
same random locations assigning them the appropriate stellar profile parameters, and magnitudes 
computed from the simulated spectra. Finally, we added the camera read noise and Poisson noise 
to the simulated data to closely mimic real sources. 

We extracted catalogs for each set of simulated quasars using the identical SExtractor param- 
eters as the unaltered images, using the i?-band image as the detection frame. We then processed 
these objects through our selection pipeline, imposing the color and morphology cuts and assign- 
ing magnitude- limits, as described in Section [2j We determine the recovered fraction of quasar 
candidates in each redshift and magnitude bin (the selection function, fsei), thereby mapping the 
completeness of our candidate selection. We list the value of fsei in column (4) of Table [5l Figure [5] 
plots contours showing the completeness as a function of redshift and magnitude for the NDWFS 
survey on the left and the DLS survey on the right; detected quasars are marked with solid circles. 

Some of our lower redshift quasars exist in areas of high incompleteness, especially in the DLS 
survey, due to the Lya line entering the R band at z ~ 4. Two of these quasars, DLS J0923-I-2952 
at z = 3.84 and DLS J1052— 0528 at z = 3.89, have unusually high Lya EWs. Since our simulated 
quasar Lya EW distribution was chosen as a Gaussian distribution with EW = 65 it 34A, there 
were no simulated spectra resembling these objects. This may point to an underlying problem with 
simulating quasars using a Gaussian distribution of EWs; on the other hand these objects may 
simply be unusual specimens. 

5.2. Determining K-corrections 

To determine the quasars' absolute luminosity, we need to perform a correction from an ob- 
served magnitude to a magnitude at a standard wavelength in the quasar rest-frame. The absolute 
magnitude at rest-frame 1450A, M1450, is most commonly used in expressions of QSO luminosity 
functions. 

Using 200 simulated spectra for each Az = 0.1 interval in range 3.5 < z < 5.2 (see Section [3]), 
we calculate the apparent magnitude m-1450 by convolving each spectrum with a boxcar filter equal 
to unity in the interval 1425A< A < 1475A and redshifted to the observed frame. The K-correction 
is defined as: 

M1450 = m - 5 log(A/10) - i^corr, (7) 

where Di is the luminosity distance in parsecs, M1450 is the absolute magnitude at rest-frame 
1450A, and m is the apparent magnitude in some fiducial filter with central wavelength A. For a 
continuum-only quasar spectrum, oc i^", the K-correction to M1450 would be equal to: 



K — 

-''-corr — 



-2.5(1 + a) log(l + z)- 2.5a log(1450A/A), 



(8) 
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where A is the effective wavelength of the filter used for the K-correction. Using simulated quasar 
spectra, the K-correction can be obtained as the difference between m and apparent magnitude 
^1^1450 in the observed frame as: 

-?^corr = m - 7711450 - 2.5 log(l + z). (9) 



We use the simulated quasar spectra to determine the expected dispersion in the value of the 
Kcorr with redshift. Figure [8] shows the value of iCcorr as well as the one sigma uncertainty interval 
determined from the various realizations of our simulated quasar spectra. The K-correction for rriR 
(black squares) sharply rises after z > 3.9, as the Lya forest progressively moves into the R band. 
At our probed redshift range, 1450A is observed at 6960 - 8990A. This suggests that although our 
selection relies on detections in the R band, our K-corrections can be based on / and z-magnitudes 
for NDWFS and DLS, respectively, minimizing the effect of emission lines on -ftTcorr- We plot the 
corrections to these bands with orange circles and red triangles. 

The blue dashed curve in . Figur e [8] shows the K-correction for the SDSS i'-band magnitude, 



determined by [Richards et al.l (|2006l ) . We also convolved our simulated quasar spectra with the 



SDSS i' filter and obtained the curve plotted w it h bla ck squares in Figure El The result is extremely 



similar to the K-correction in [Richards et al.l (|2006l ) below z = 4.8, after which the inclusion of 



different descriptions for Lya emission and the Lya forest simulations in our spectra create a 
progressive offset. The / and z-band corrections are both consistent with the SDSS corrections out 
to z = 4.3 after which they diverge slightly, but never as much as for the i?-band. The dispersions 
in the K-corrections for the SDSS nii and our surveys' magnitudes are equal, having a mean value 
of CTKcorr = 0.08 in both cases. 

Because our quasars are at z ~ 4 and all our spectra cover 1450A in the quasars' rest frame, we 
can also measure M1450 directly from the spectra. To do this, we compute the magnitude from the 
median flux, 1450, between res t-frame A = 1425A and A = 1475A, mi45o = —2.41 — 2.5 Iog(/A,i45o) 



(Eqn. (2) of lBlanton et al.ll2003l ). We also determine the spectrophotometric magnitude, -Rgpecphot, 
of each object by convolving its spectrum with the i?-band filter profile and scaling it to the image- 
based i?-band magnitude. We determine the exact K-correction for each spectrum, independent of 
any modeling, where 

Kcorr = -Rspocphot " "T-MSO " 2.5 log(l + z), (10) 

and apply the K-correction as shown in Equation [71 

This in principle should give us a more accurate measure of M1450 over the simulated K- 
corrections. However, the wavelength range of our LRIS spectra spans 5750-8300A, while the R 
filter profile has 64% transmission at 5750A and drops to 0% below 5000A (see Figured]). We 
could, therefore, be underestimating the flux transmitted through the i?-band, especially for the 
lower-redshift objects in our sample. 

Figure [9l plots the comparison between M1450 derived using our simulated spectra (plotted on 
the vertical axis) and directly from the object spectra (plotted on the horizontal axis). The triangles 
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are quasars from NDWFS and the asterisks are quasars from DLS. The dotted hne represents a one- 
to-one correlation, and while the scatter about this line can be significant (up to ~ 0.3 magnitude), 
the best fit line to these points (plotted as a solid black line) has a slope of 1.02 it 0.09 with an 
intercept of 0.19 ± 2.06, suggesting an offset of ~ 0.3 magnitudes. The M1450 measured directly 
from the spectra tend to be fainter than M1450 estimated from the simulations, which may (at least 
in part) be because of the incomplete spectroscopic coverage of the i?-band. 

Table [5] lists the M1450 magnitudes computed from applying the K-corrections from our sim- 
ulated spectra and from the object spectra directly in columns (6) and (9), respectively. In the 
following sections, we use both of these values to compute and analyze the QLF. 



5.3. The z ~ 4 Quasar Luminosity Function 



With the selection function and absolute m agnitudes in ha nd, we are able to calculate the QLF 



at 2: ~ 4. The 1/Va (also known as the WV^c^: ISchmidtill968l ) method is the most straightforward 
way to compute the volume density of quasars as a function of luminosity, ^{Mi4^q,z = 4). The 
available volume, Va, is the comoving volume within which a quasar with redshift z and magnitude 
R can be found. We determine this volume in pencil-beam units of Mpc^ arcmin"^ and multiply 
by the area surveyed, A, of the NDWFS and DLS, 1.71 deg^ and 2.05 deg^, respectively, correcting 
by our selection function and spectroscopic incompleteness (i.e., the fraction of candidates without 
spectra) . 



5.3.1. The R<23 Quasars 

For candidates with R < 23, our spectroscopic completeness is 0.75 and 0.71 in NDWFS and 
DLS, respectively. Beyond this, the completeness drops significantly to 0.33 (2/6) and 0.11 (1/8) 
for 23 < R < 23.5 (and further, stih, for R > 23.5 to 0.04 and 0.02) in the NDWFS and DLS, 
respectively. This effectively places our survey limit at i? = 23 and we compute here the luminosity 
function for only these objects. We divide our sample into one-magnitude- wide bins, which span 
—27 < M1450 < —20, for each survey separately, and together, as a combined sample. 

While our survey sought to find quasars with 3.8 < z < 5.2 and brighter than R = 23.0, 
corresponding to 4291 Mpc^ arcmin"^, this is not always the volume available to the quasars that 
we find. For example, a quasar with R = 22.9 at z = 4.1 has M1450 = —23.24. The limiting redshift 
for such a quasar would be, ziim = 4.18, beyond which it would be too faint to make our magnitude 
limit. Therefore, Va for this quasar is the comoving volume between z = 3.8 to z^^, which, in this 
example, is 1427 Mpc^ arcmin"^. We calculate ziim and Va for each of our quasars separately, and 
list them in Table [5l 

The selection function, /sei, which we described in §5.1^ is the probability that a quasar with 
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magnitude R and redshift z would be selected as a candidate. For each quasar, we scale Va by 
the selection function and spectroscopic completeness, /spec (Column (7) of Table [5] referring to 
the fraction of candidates with spectroscopic followup), in each luminosity bin. We compute the 
volume density of quasars for each M1450 bin for NDWFS and DLS separately using the following 
formalism: 

N 

$(Mi450, Z = A) = ^(/sel(i) X /spec(i) X Va{i) X A)"', (11) 
i 

where N is the number of quasars in the luminosity bin. We list the values of ^ for NDWFS and 
DLS as well as their average, in Table El Figure fTOl plots <I>(Mi45o,z = 4), using M1450 from K- 
corrected z and / magnitudes on the left and directly from the quasar spectra on the right. We plot 
the volume densities of NDWFS (squares) and DLS (triangles) quasars as well as for the combined 
sample (filled circles). We fit a single-power-law, <1?(M) = cj)*io~''-^(^~^)*^, to these measurements. 
We find that the shape of the QLF depends strongly on how we compute M1450. When we compute 
M1450 using the K-corrected z and / magnitudes, the volume density resembles a power law (with 
a slope of /3 = —1.59 it 0.22) out to the faintest bin, which suffers from the flux-limit bias of our 
survey (which we discuss in the next subsection) . When we compute M1450 directly from the quasar 
spectra the objects populate fewer bins. While the shape of this QLF is scattered and does not 
resemble a power-law, the space densities derived from the individual surveys are more consistent 
with each other. Table [7| lists our fitted value for (3 from both datasets. 



5.3.2. The R > 23 Quasars 

Of the 148 candidates that we selected in ^3.21 117 have R > 23, 59 in the DLS survey and 
63 in the NDWFS survey. We obtained spectra for six of these objects; two NDWFS sources were 
Lyman break galaxies (LBGs) at z ~ 3.8, three were quasars in our desired redshift range (two in 
the DLS and one in the NDWFS) and one NDWFS source was a Type-II AGN at z = 0.21. We 
wish to make use of the three i? > 23 quasars to enhance our measurement of the QLF. 

We compute the volume density for these quasars according to Equation [TTJ However, we 
determine the limiting redshift, z\[^, and Va assuming a survey magnitude limit of R = 23.5 
and R = 24, when appropriate. The largest source of uncertainty in this method arrives from 
spectroscopic completeness, /spec, which we use to scale our available volume Va and which suffers 
from small number statistics. For example, in the DLS sample, we obtained one spectrum (of a 
quasar) out of nine candidates with 23 < i? < 23.5, rendering /spec = 5 = 0.11 and the Poisson 
error on this is 0.12. 

These values are listed in Table [5] and we compute the QLF with these three additional quasars 
using both estimates of We plot the QLF including all of our quasars in Figure [TTJ Again, 

we fit a single power-law to these data and find that regardless of how we compute M1450, the QLF 
resembles a power-law. The slopes are marginally consistent, and steeper than the fit using only 
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quasars with i? < 23. This is not surprising, as the volume densities for the faintest quasars are 
potentiahy significantly inflated because of the low spectroscopic completeness, which appears in 
the denominator of Equation llli Nevertheless, including these data appear to remedy the divergent 
volume densities computed for NDWFS and DLS in the left-hand panel of Figure [TO] as well as the 
scattered shape of the QLF in the right-hand panel of Figure [TOl 

We list the volume density for our QLF including all 23 quasars for both methods of computing 
-^1450 in Table [6] and the best-fit slope in Table [71 



5.4. Comparison with other surveys 



We compare our binned QLF data points for 3.8 < z < 5.2 (median z = 4.15) with the results 
from other surveys at similar redshifts in Figure [I2j Our results for the binned QLF using all the 
quasars in our survey, including those with R > 23, are plotted with red circles. The left-hand 
panel shows our QLF where M1450 has been computed from K-corrections to the z and / band 
photometry and the right-hand panel shows the QLF where M1450 was obtained directly from the 
spectra. We plot with a dotted line the best-fit power-law for the faint end, shown in Figure [TTl 

To extend our Q LF into the bright end, the most suitable sample comes from the SDSS. 
Richards et al.l (l2006l ) determined the QLF for SDSS quasars from the Third Data Release (DR3). 
We utilize their results for the z = 4.25 bin and plot their binned LF with asterisks and their 
single-power-law model fit with a dash-dot-dot-dot line in Figure [121 Both the po i nts an d the curve 
are evolved to 2; = 4.15 using the evolution model determined by 



Eichards et al.l pnOffl . 



Fontanot et al.l (120071 ) present a re-analysis of the SDSS-DR3 quasars with M1450 < —26 



combined with eleven quasars from the Great Observatories Origins Deep Survey (GOODS) with 
M1450 > —23.5 in the range 3.5 < z < 5.2. We plot their QLF at z = 3.75 and z = 4.6 with open 
triangles and squares in Figure US We evolve both these data sets to z = 4.15 using the sarn e 
evolution model as above. The pure density evolution (PDF) model fit from lFontanot et al.l (|2007l ). 
also evolved to z = 4.15, is plotted with a dash-dot line. 



While [Richards et al.l (|2006l ) and [Fontanot et al.[ (|2007[ ) rely on the same SDSS-DR3 quasars, 
their derived volume densities are highly discrepant. The difference between the LFs s tems from 
their d ifferent estimates for the completeness of the SDSS sample as a function of redshift. [Fontanot et al, 
([20071 ) estimate the completeness for the SDSS subsample use d in their calcul a tion b y building a 
spectral library based on the quasar template spectrum from [Cristiani &: Viol (|l990l ). Instead of 
modifying t he composite using a Gaussian distribution of power-law continua (as is done in this 
work and in lRichards et allbood ). they empirically determine a distribution of continua blueward of 
Lya from 215 high quality quasar spectra with 2.2 < z < 2.25 in the SDSS-DR3 quasar sample. As 



a consequence, the distribution of power-law slopes in the [Fontanot et al.[ ([2007 ) library is steeper 



with 02 



-0.7 ± 0.3, compared with 02 



-0.5 ± 0.3 in the 



Richards et al 



(|2006l ) treatment. 



The IGM absorption is then added in a manner similar to ours. They use this spectral library to 
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estimate the completeness of the SDSS color selection criteria ([Richards et al.ll2002l ) and find that 
the completeness drops to 50% (and as low as 30% at z = 4.4) for 3. 6 < z < 4.4. At 4.5 < z < 4.9 
the completeness rises to above 90% and drops off again after (see iFontanot et al.l 120071 . Figure 
A.4.II 



Using a method similar to ours to create their spectral library, iRichards et al.l (120061 ) calcu- 
late a completeness o f nearly 90% for z > 3.8. While our method m ore closely matches th at of 



Richards et al.l (|2006l ) , our values seem to have better continuity with iFontanot et al.l (j2007l ) , es- 
pecially for the QLF whose M1450 was derived from K-corrected z and / band magnitudes (the 



left-hand panel of Figure [T2]) . 

For comparison, we overplot in Figure [12] QLF model fits for several surveys at or near our 
survey's redshift. The left-hand panel shows general agreement between our QLF based on M1450 
computed froni simul ated quasars (in all but the fa intest bins) and the SD SS+GOODS analysis of 



Fontanot et al.l ([2007[ ) as well as the VVDS survey ([Bongiorno et al.[[2007l ). The latter spectroscop- 
ically identified 130 broad-line AGN with z = — 5 over ~ 1.7 deg2 in two fields. The z = 3.0 - 4.0 
bin from their QLF contains 17 QSOs with a luminosity range of —26 < M1450 ^ —22 and is 
well-suited for comparison with our data. We plot their model fit with a solid black line in Figure 

m 

There is a large discrepancy between our m easurements and predictions for the value of the 
quasar density at z ~ 4. The COMBO-17 survey (jWolf et ahlbood ) covered 0.78 deg^ to < 24 and 
found 11 QSOs in the redshift range 3.6 < z < 4.2, spanning the luminosity range —27 < M1450 < 
—24 and 4 QSOs in the redshift range 4.3 < z < 4.8 in the luminosity range of —27 < M1450 < —26. 
While their results are in agreement with SDSS at the bright end, their QLF modeled as Pure 
Density Evolution (PDF, dashed line in Figure [T2|) significantly underpredicts the number of quasars 
at the faint end. 



5.5. Model Fitting to the Combined QLF 



The shape of the QLF is typically parametrized by a standard double power-law form: 



$(M, z) 



$(M*) 

]_Q0.4(a+l){M-M*) _j_ ]_Q0.4(/3+l){Af-M*) 



(12) 



We follow the formalism of lSandage et al.l ([l979l ) and lMarchesini et al.l ([20071 ) to determine the best 
fit parameter values for the QLF using the maximum likelihood (ML) method. We combine the 
SDSS quasars at 2: ~ 4 with our 23 quasars to sample as broad a magnitude range as possible when 
fitting the double-power-law in equation [T2l We try dif ferent combinat i ons o f our QLF, based on 
the two ways of calculating M1450 , with the QLF from [Richards et al.l ([20061 ) and IFontanot et al 



^This is discussed in detail in appendix (A. 3) of IFontanot et al.l ( 20071 ) 
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(j2007l ) at the bright end. As is clear from Figure [121 the bright-end QLF from lFontanot et al.l (|2007l ) 
has better continuity with our QLF b ased on Mi 450 c ompu ted from simulated quasars. There is 
a sharp jump between our values and [Richards et al.l (|2006l ). We therefore do not attempt to fit 
these points in this case. We do not attempt to model the evolution of the QLF, since our survey 
is focused on a "snapshot" of the quasar population at a narrow redshift range. 

We initially attempted to maximize the likelihood A = n£iP« with respect to all three 
parameters, a, (3 and M* . We created a gridded cube of a, /3 and M* values, spanning a, /3 = [—4, 0] 
and M* = [—27,-20], and computed A at each point, first using a "coarse" grid (Aa,/3 = 0.1, 
AM* = 0.2) with the intention of refining the values near the maximum. However, due to the 
small number of quasars in our sample, poorly constraining the faint end, no clear maximum could 
be found over these reasonable ranges of a, (5 and M*. 

Since the ML method fails to constrain the double-power-law parameters, we fit Equation [T2l 
to the binned QLF using a minimization and allowing all four parameters, <!>*, M*, a, and (3, 
to vary. The results of our fit are listed in Table [8] for the three combinations of data sets that we 
fit. We plot with a red line in Figure [12] the best-fit double-power-law parameters from dataset (1) 
in the left-hand panel and dataset (3) on the right-hand panel. 

When our QLF is combined with SDSS data on the bright end, the best-fit parameters are 
more strongly influenced by the choice of bright-end dataset than our method for estimating M1450. 
This is because the error bars on the bright-end bins are much smaller compared to the faint-end 
and strongly constrain the fi t. Furthermore , despite the different shapes of the QLF in our two 
derivations, the fit using the iFontanot et al.l (|2007l ) are nearly identical and largely unaffected by 
the four bins made up of GOODS quasars on the faint end; we excluded these bins and obtained 
nearly the same values for the double-power-law parameters. 



When fitting to our data combined with [Richards et al.l ([20061 ) . we compute, effectively, a 
single power law fit, with a slope a ~ /3 = —2.3 ± 0.2. This result is provocative, both because it 
preserves the slopes measured for each data set individually (q ~ 2.1 for z = 4.25 in Figure 21 of 
Richards et al.|[2006[ . while we find (3 = —2.46 in the right-hand panel of Figure [TT]) . and because 
this suggests that the number of faint AGN rises dramatically out to very faint magnitudes, and 
that there is no observed "knee" to the QLF at z ~ 4. This degeneracy between a, (5 and M* 
explains why the ML method for the double-power-law fit did not converge. Therefore, we fit a 
single power-law, $(M) = <^*\Q~'^-^'yP~^)^^ ^ to this data set using the ML method. In this case 
there is one free parameter, /?, which we allow to vary between —1 and —4 in steps of A/3 = 0.01. 
The likelihood is maximized at (3 = —2.67 and the 68% confidence limits are (-2.80,-2.55), the 
90% confidence Umits are (-2.88, -2.48) and the 99% confidence limits are (-3.00, -2.37). While 
this result stands in conflict with other measur ements of the QLF to faint magni tudes which see 



a flattening of the faint end with redshift (e.g., [Wolf et al 
plausibility and possible interpretation in ^ 



2003 



Hunt et al.[[2004[ ). we discuss its 
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6. Estimating the Contribution of Quasars to the UV Radiation Field at z ~ 4 



We use the best-fit parameters from our minimization to all the binned data to integrate 
Equation [T2l and compute the emissivity of quasars at 1450A, €1450- Our computed UV radiation 
depends strongly on the best-fit parameters (listed in Table [8]) for the QLF, especially the faint-end 
slope, [3. 

For dataset (1), wh ere we combined our b inned QLF using K-corrected z and /-band photom- 
etry with the data frorn iFontanot et alj (|2007l ). we compute £1450 = 9 x 10^^ erg s^^ Hz^^ Mpc^'^. 
Following iMadau et al.l (jl999l ). we use our parametrized QLF together with the quasar SED of 



Elvis et al. 



(|1994| ) to compute the photoionization rate from QSOs at this redshift, Nqso = 4 x 10^^ 
s~^ (integrating J ^{L)LdL down to M1450 = —20). This is almost twice the value needed to ionize 
the intergalactic medium (IGM) at z = 4.15, A'igm = 2.4 x 10^^ s~^ and stands in c ontrast to pre- 



vious statements on the abil i ty of quasars to ionize th e universe at higher redshifts ([Haiman et al 



2001 



Wvithe Loeb 



2003a 



Shankar &: Mathuijl2007l ) which suggested that AGN do not produce 



a sufficient number of ionizing photons and starforming galaxies must play a larger role. 

For datas et (3), where we comb ined our binned QLF using M1450 from the object spectra with 
the data from [Richards et al.l ( 20061 ) and effectively fit a single power law, we find £1450 = 7 x 10^^ 
erg Hz~^ Mpc~^ and Nqso = 3 x 10^^ photons s~^ Mpc^^, which also produces sufficient 
photons to ionize the IGM. This value is strongly dependent on the faint limit of integration. 
Because the total UV-luminosity density has to be larger from the contribution of quasars fainter 
than our integration limit, our computation of Nqso is a lower limit of the photoionization rate. 



7. Discussion 

Our survey has found that the comoving volume density of quasars continues to rise as a steep 
power law to low luminosities, four magnitudes fainter than previous measurements. This is true 
regardless of how M1450 is computed. This result conflicts with predictions for the evolution of 
the shape of the QLF with redshift, based on the observed evolution of the QLF at lower z. In 
addition, this is surprising because objects at the faint end of a flux-limited survey, with their 
associated large photometric errors, tend to be undercounted. While our results are sensitive to 
the fluctuations of small number statistics (our faintest bin contains between one and two objects, 
depending on how M1450 is computed) it behooves us to come up with a reasonable explanation or 
rule out any systematics that may lead to an excess of quasars counted in error. 

In this section, we consider possible explanations for overcounting our quasars at the faintest 
bins. These are: (1) we are including galaxies in our sample, (2) our selection function overcorrects 
$(M) (equation II ip leading to an overestimate of the QLF, (3) clustering of quasars in our chosen 
flelds due to large scale structure. We examine each of these possibilities below. 

We rule out possibility (1), that we have counted galaxy interlopers as quasars, by an examina- 
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tion of FigureO All of the spectra except for our lowest-redshift object, NDWFS J142713. 2+322842 
at z = 3.74, show strong, broad Lyman-a emission with the classic asymmetric profile from 
absorption of the blue wing of the emission line attributed to high-redshift quasars. NDWFS 
J142713. 2+322842 is also obviously a quasar, as its C IV line is extremely broad, with FWHM 
~ 11,000 km s"^ 

To examine possibility (2), that our number counts are too high by an underestimate of the 
selection function, we plot the value of the selection function for each quasar versus its absolute 
magnitude, M1450 in Figure [13l As we already noted in Section [5TH fsei is ~ 80% for most of the 
quasars. At faint luminosities, f^d drops to ~ 0.4, which means that <I>(Mi45o,2 = 4) in equation 
[TTjis corrected for incompleteness by a factor of ~ 2.5. This is insufficient to account for the steep 
rise in $(M) that we see, since a flattening of the faint-end slope of the QLF would predict fewer 
quasars by a factor > 50. 

The third possibility is that clustering due to cosmic variance is enhancing the number of 
quasars that we find, which can have a significant impact in a small sample. We examined the 
spatial and redshift distribution of our quasars for each field and found that of the three faintest 
and lowest luminosity quasars (which dominate the last two bins in our QLF) are isolated, with 
the nearest quasar > 0.5 deg away. This corresponds to a separation on the sky of ~12 Mpc at 
z ~ 4.15. 

DLS J105346. 1—052859, the faintest quasar in our sample with R = 23.83, is 4f4 away from 
DLS J105348. 8— 053319, a much brighter quasar with R = 20.94. While these objects are the 
closest pair of quasars in our sample, their redshift difference (the faint quasar is at z = 4.02 and 
the brighter quasar is at z = 4.20, implying that their orthogonal separation is ~ 2000 Mpc in our 
stated cosmology) rules that their proximity is merely a projection effect. 

There ar e a few qua s ar sa mples that, while small, are consistent with steeper QLF at fainter 
luminosities. ICool et al.l (|2006l ) found three z > 5 quasars in NDWFS Bootes using mid- infrared 
color selectio n, which is mearit to be insensitive to dust reddening and avoids confusion with the 
stellar locus (jStern et al.ll2005l ). None of these quasars appear in our candidate list from NDWFS. 
This is because our color selection be comes highl y incom plete at z > 5 and R > 23. In addition, 
the optical colors of the quasars from ICool et al.l (|2006l ) fail to meet our color criteria. Only one 
object, J142937. 9+330416, has Bw — R > 2.2, and none of them met the criteria set forth in eqn. 
(2). While their objects are somewhat more luminous than our z ~ 4 quasars, with M1450 = —26 
to —24.5, they are 1.5 magnitudes fainter than the SDSS quasars at these redshifts. The space 
density derived from the discovery of these three quasars is consistent with a luminosity function 
that has a steep power-law index, ^{L) oc discovery of a r adio-loud, z = 6.12 quasar in 

NDWFS (the probability of which is < 1%) by iMcGreer et al.l (j2006l ) also suggests that the space 
density of faint, high-redshift quasars may be higher than previously thought. 
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8. Conclusion 



Using 23 quasars at z ~ 4 discovered in our survey of deep optical imaging data from the DLS 
and NDWFS surveys, we have measured the faint end of the QLF. Depending on how we compute 
Mi45o, our directly fit faint-end slope ranges between /? = —1.98 =fc 0.23 and /3 = —2.46 it —0.20, 



in both cases steeper than the faint-end slope measured by iFontanot et al.l (|2007l ) who find /3 = 
—1.71 lb 0.41. When we combine our QLF at the faint end with with the SDSS-based QLF at the 
bright end from Richards et al. ( 20od ) we conclude that the shape of the z ~ 4 QLF is best fit by 
a single power law, with a slope of ~ —2.7 it 0.1 (via the ML method). A QLF with this shape is 
able to produce enough UV photons to ionize the IGM at this redshift, which is unexpected. 

This result is pro vocative and presents a challenge to interpretat ions of the shape of the QLF 



rp 
1 



and its evolution (e.g. iHopkins et al.ll2006al : IShankar &: Mathurj|2007l ) and has two immediate cos- 
mological implications: (1) the models of the faint AGN evolution at the epochs ~ 1 Gyr after the 
end of the reionization, and possibly all the way into the reionization era, may need to be revisited; 
(2) AGN were a more significant contributor to the metagalactic ionizing UV flux at these epochs, 
affecting the evolution of the IGM. We caution, however, that this result is currently reliant on 
only a few quasars at the faintest end of our survey. Additional spectra of our i? > 23 candidates 
will provide a more robust measurement in this crucial luminosity regime and better constraints 
for theoretical models. 
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Table 1. NDWFS Fields Surveyed 



Field 


R.A. 


Dec 


Size 


FWRMbw 


FWHMfl 


FWHMj 


EXPbw 


EXPr 


EXP I 




(J2000) 


(J2000) 


(0 


n 


(") 


(") 


(sec) 


(sec) 


(sec) 


NDWFS J 1426+3236 


14:26:03.74 


+32:36:31.72 


36.8 X 38.0 


0.93 


0.94 


0.96 


8400 


6000 


12000 


NDWFS J1431+3236 


14:31:36.14 


+32:36:46.29 


36.8 X 38.0 


0.80 


0.91 


0.77 


8400 


6000 


12000 


NDWFS J1434+3421 


14:34:30.79 


+34:21:54.18 


35.8 X 37.4 


1.05 


0.98 


0.79 


8400 


4200 


10200 


NDWFS J1437+3347 


14:37:16.32 


+33:47:01.72 


36.9 X 38.0 


0.89 


0.87 


0.88 


8400 


6000 


12000 


NDWFS J1437+3457 


14:37:24.59 


+34:57:02.13 


35.5 X 38.0 


0.86 


1.07 


1.16 


8400 


10800 


11400 



Table 2. DLS Fields Surveyed 



Field 


R.A. 


Dec 


Size 


FWHMb 


FWHMy 


FWHMij 


FWHM^ 


EXPr 


EXP bVz 




(J2000) 


(J2000) 


(') 


(") 


(") 


(") 


(") 


(sec) 


(sec) 


F2p21 


09:22:37.1 


+30:00:00 


35.1 X 35.1 


1.34 


1.08 


0.86 


0.88 


18000 


12000 


F2p31 


09:22:37.1 


+29:20:00 


35.1 X 35.1 


1.26 


1.05 


0.78 


1.50 


18000 


12000 


F2p33 


09:16:27.7 


+29:20:00 


35.1 X 35.1 


1.24 


1.02 


0.86 


1.28 


18000 


12000 


F4p23 


10:49:19.4 


-05:00:00 


35.1 X 35.1 


1.09 


1.10 


0.89 


1.12 


18000 


12000 


F4p31 


10:54:40.8 


-05:40:00 


35.1 X 35.1 


1.20 


0.99 


0.87 


1.11 


18000 


12000 


F4p32 


10:52:00.0 


-05:40:00 


35.1 X 35.1 


1.17 


0.93 


0.89 


1.24 


18000 


12000 



Table 3. NDWFS Candidates 



R.A. 


Dec 


Field 


Bw 


R 


/ 


— a 


R-I 


FWHMfl 


Redshift 


Class. 


(J2000) 


(J2000) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(arcsec) 






14 25 22.71 


+32 48 27.0 


1426p3236 


23.35+0.02 


21.10+0.01 


20.96+0.01 


2.25 


0.14 


1.08 


3.81 


QSO 


14 30 17.78 


+32 20 03.4 


1431p3236 


>27.90 


22.23+0.01 


21.43+0.01 


5.68 


0.80 


1.03 


4.75 


QSO 


14 30 35.22 


+32 44 26.3 


1431p3236 


24.08+0.03 


21.49+0.01 


21.50+0.01 


2.59 


-0.01 


0.99 


0.21 


AGN II 


14 31 01.94 


+32 51 55.5 


1431p3236 


25.52+0.11 


21.77+0.01 


21.17+0.01 


3.75 


0.60 


0.97 


3.94 


QSO 


14 31 52.32 


+32 54 50.2 


1431p3236 


23.94+0.04 


21.42+0.01 


21.12+0.01 


2.51 


0..30 


1.03 






14 32 04.74 


+32 54 05.8 


1431p3236 


24.61+0.05 


22.31+0.01 


22.10+0.01 


2.30 


0.21 


0.97 






14 33 24.54 


+34 08 41.2 


1434p3421 


26.13+0.19 


22.62+0.02 


22.61+0.02 


3.51 


0.01 


0.92 


3.88 


QSO 


14 33 31.15 


+34 32 48.3 


1434p3421 


24.24+0.04 


20.52+0.00 


20.53+0.00 


3.71 


0.00 


0.92 


4.15 


QSO 


14 36 26.69 


+34 49 50.1 


1437p3457 


25.38+0.10 


22.93+0.02 


22.73+0.03 


2.45 


0.20 


1.13 






14 36 42.86 


+35 09 23.8 


1437p3457 


24.66+0.05 


22.00+0.01 


21.95+0.02 


2.66 


0.05 


1.17 


3.90 


QSO 


14 36 58.34 


+33 36 32.0 


1437p3347 


23.72+0.02 


20.43+0.00 


20.31+0.00 


3.29 


0.12 


0.93 


4.02 


QSO 


14 37 32.67 


+33 55 22.0 


1437p3347 


26.61+0.35 


22.88+0.02 


22.57+0.03 


3.74 


0.30 


0.96 


4.22 


QSO 


14 37 34.26 


+34 53 32.9 


1437p3457 


25.39+0.10 


22.65+0.02 


22.33+0.02 


2.74 


0.33 


1.14 




? 


14 37 50.50 


+34 59 52.9 


1437p3457 


25.49+0.11 


22.94+0.02 


22.63+0.03 


2.55 


0.31 


1.27 






14 38 13.85 


+35 02 36.4 


1437p3457 


25.75+0.13 


22.85+0.02 


22.67+0.03 


2.90 


0.18 


1.26 


4.25 


QSO 


14 38 39.68 


+35 12 45.9 


1437p3457 


24.69+0.05 


21.10+0.00 


20.52+0.01 


3.59 


0.58 


1.12 


4.63 


QSO 



Table 4. DLS Candidates 



R.A. 


Dec 






Field 


B 


V 


R 


z 


B- R 


B-V 


R- z 


FWHMfl 


Redshift 


Class. 


(J2000) 


(J2000) 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 








09 15 27.53 


+29 17 


50. 


,4 


F2p33 


24.73+0.07 


21.76+0.01 


20.72+0.00 


20.16+0.01 


4.02 


2.97 


0.55 


0.89 


4.34 


QSO 


09 21 24.62 


+29 53 


30. 


,8 


F2p21 


24.91+0.09 


23.42+0.03 


22.57+0.01 


23.28+0.13 


2.34 


1.49 


-0.71 


0.87 






09 21 51.96 


+29 24 


57. 


,2 


F2p31 


25.41+0.14 


23.44+0.03 


22.46+0.01 


21.86+0.04 


2.95 


1.97 


0.60 


0.81 


4.32 


QSO 


09 22 36.49 


+30 10 


10, 


,6 


F2p21 


25.53+0.16 


24.23+0.07 


22.90+0.01 


22.94+0.10 


2.63 


1.30 


—0.04 


0.88 


3.98 


QSO 


09 23 27.22 


+29 52 


51, 


,7 


F2p21 


23.82+0.04 


22.28+0.01 


21.58+0.01 


21.66+0.03 


2.23 


1.54 


-0.08 


0.97 


3.84 


QSO 


09 23 36.82 


+30 09 


49, 


,9 


F2p21 


>27.36 


24.28+0.07 


21.94+0.01 


20.83+0.02 


5.42 


3.08 


1.11 


0.90 


5.06 


QSO 


10 51 19.41 


-05 55 


25, 


,7 


F4p32 


24.17+0.05 


22.74+0.02 


21.94+0.01 


22.14+0.04 


2.23 


1.43 


-0.20 


0.93 






10 51 54.74 


-05 46 


26. 


,8 


F4p32 


24.87+0.09 


22.97+0.02 


22.01+0.01 


21.70+0.03 


2.86 


1.90 


0.31 


1.03 






10 52 19.78 


-05 28 


18, 


,2 


F4p32 


25.35+0.14 


23.66+0.03 


22.92+0.01 


22.92+0.08 


2.43 


1.70 


0.00 


0.93 


3.89 


QSO 


10 52 27.95 


-05 42 


34, 


,7 


F4p32 


22.84+0.02 


20.64+0.00 


19.68+0.00 


19.13+0.00 


3.16 


2.19 


0.55 


0.88 


3.89 


QSO 


10 53 48.89 


-05 33 


19, 


,4 


F4p31 


23.82+0.03 


21.65+0.01 


20.94+0.00 


20.54+0.01 


2.88 


2.17 


0.40 


0.99 


4.20 


QSO 


10 55 07.12 


-05 30 


14, 


,9 


F4p31 


25.61+0.15 


22.37+0.01 


21.02+0.00 


20.35+0.01 


4.59 


3.23 


0.67 


0.93 


4.40 


QSO 


10 55 23.03 


-05 48 


50, 


,7 


F4p31 


26.27+0.27 


23.34+0.03 


22.35+0.01 


21.82+0.03 


3.92 


2.92 


0.52 


0.90 


4.12 


QSO 


10 55 44.41 


-05 31 


55, 


,9 


F4p31 


24.44+0.05 


23.11+0.02 


22.06+0.01 


22.71+0.06 


2.38 


1.34 


-0.65 


0.96 







Table 5. z ~ 4 Quasar Sample 



Simulated, z, I Spectrum 

Name R z fsei /spec M1450 ^um'' Va" M1450 aiim°' K'' Area 

(mag) (mag) (Mpc^ arcmin"^) (mag) (Mpc^ arcmin"^) (deg^) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



NDWFS 



NDWFS J142522. 71+324827.0 


21.10 


3.81 


0.723 


0.75 


-24.74 


3.81 


4024.4 


-23, 


.79 


3.81 


2588.5 


1.71 


NDWFS J142713. 19+322840.8 


23.70 


3.74 


0.387 


0.036 


-22.17 


3.74"= 


1248.8"= 


-21, 


.46 


3.74"= 


278.5"= 


1.71 


NDWFS J143017.78+322003.4 


22.23 


4.75 


0.655 


0.75 


-24.80 


4.75 


4104.7 


-23, 


.89 


4.75 


2786.6 


1.71 


NDWFS J143101.94+325155.5 


21.77 


3.94 


0.806 


0.75 


-24.60 


3.94 


3862.9 


-24, 


.90 


3.94 


4238.0 


1.71 


NDWFS J143324.54+340841.2 


22.62 


3.88 


0.708 


0.75 


-23.13 


3.88 


1159.2 


-22 


.63 


3.88 


278.5 


1.71 


NDWFS J143331. 15+343248.3 


20.52 


4.15 


0.842 


0.75 


-25.37 


4.15 


4291.1 


-25, 


.54 


4.15 


4291.1 


1.71 


NDWFS J143642.86+350923.8 


22.00 


3.90 


0.782 


0.75 


-23.80 


3.90 


2616.9 


-23, 


.29 


3.90 


1516.1 


1.71 


NDWFS J143658.34+333632.0 


20.43 


4.02 


0.839 


0.75 


-25.51 


4.02 


4291.1 


-25 


.26 


4.02 


4291.1 


1.71 


NDWFS J143732.67+335522.0 


22.88 


4.22 


0.724 


0.75 


-23.36 


4.22 


1692.9 


-23 


.01 


4.22 


858.4 


1.71 


NDWFS J143813.85+350236.4 


22.85 


4.25 


0.726 


0.75 


-23.28 


4.25 


1486.5 


-22, 


.95 


4.25 


737.2 


1.71 


NDWFS J143839.68+351245.9 


21.10 


4.63 


0.805 


0.75 


-25.65 


4.63 


4291.1 


-25, 


.15 


4.63 


4291.1 


1.71 


DLS 


DLS J091527.53+291750.4 


20.72 


4.34 


0.909 


0.714 


-25.84 


4.34 


4291.1 


-25, 


.60 


4.34 


4291.1 


2.05 


DLS J092151. 96+292457.2 


22.46 


4.32 


0.805 


0.714 


-24.13 


4.32 


3178.3 


-23 


.75 


4.32 


2531.6 


2.05 


DLS J092236.49+301010.6 


22.90 


3.98 


0.497 


0.714 


-22.86 


3.98 


462.8 


-22, 


.11 


3.98 


278.5 


2.05 


DLS J092327.22+295251.7 


21.58 


3.84 


0.090 


0.714 


-24.05 


3.84 


3039.1 


-23, 


.11 


3.84 


1129.2 


2.05 


DLS J092336.82+300949.9 


21.94 


5.06 


0.739 


0.714 


-25.54 


5.06 


4291.1 


-25, 


.22 


5.06 


4291.1 


2.05 


DLS J105011.52-044253.9 


23.24 


4.27 


0.737 


0.125 


-23.07 


4.27'' 


2158.9'' 


-22, 


.70 


4.27'' 


1338.2'' 


2.05 


DLS J105219.78-052818.2 


22.92 


3.89 


0.228 


0.714 


-22.82 


3.89 


340.0 


-22 


.51 


3.89 


278.5 


2.05 


DLS J105227.95-054234.7 


19.68 


3.89 


0.302 


0.714 


-26.61 


3.89 


4291.1 


-25, 


.52 


3.89 


4291.1 


2.05 


DLS J105346.14-052859.5 


23.83 


4.02 


0.371 


0.021 


-20.40 


4.02"^ 


278.5"= 


-21, 


.74 


4.02"= 


278.5"= 


2.05 


DLS J105348.89-053319.4 


20.94 


4.20 


0.924 


0.714 


-25.38 


4.20 


4291.1 


-24 


.78 


4.20 


4078.0 


2.05 


DLS J105507.12-053014.9 


21.02 


4.40 


0.891 


0.714 


-25.68 


4.40 


4291.1 


-25, 


.18 


4.40 


4291.1 


2.05 



Table 5 — Continued 



Simulated, I Spectrum 
Name R z fsei /spec M1450 znm'^ V^^ M1450 «iim^ Va* Area 
(mag) (mag) (Mpc^ arcmin"^) (mag) (Mpc'' arcmin"^) (deg'^) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (H) (12) 

DLS J105523.03-054850.7 22.35 4.12 0.779 0.714 -24.06 4.12 3067.0 -24.40 4.12 3564.0 2.05 



*These values are computed to a survey magnitude limit oi R = 23, unless otherwise specified. 
''These values were computed to a survey magnitude limit oi R = 23.5. 
■^These values were computed to a survey magnitude limit oi R = 24. 



Table 6. Luminosity Function 





NDWFS 




DLS 






All Quasars 




-^1450 bin center 








-^Vqso 


< Ml450 > 


$ 


Nqso 


(mag) 


(10"'^ Mpc"^ mag" 


1) (10 


Mpc"^ mag"^) 




(mag) 


(lO"'^ Mpc-3 mag-^) 










Mi4B0 from simulated quasars and 


0, / magnitudes 






-26.5 






1.46±1.46 


1 


-26.61 


0.80±0.80 


1 


-25.5 


1.83±1.06 


3 


2.06±1.03 


4 


-25.58 


1.95±0.74 


7 


-24.5 


2.25±1.30 


3 


8.48±7.03 


3 


-24.44 


5.65±3.88 


6 


-23.5 


7.47±3.90 


4 


6.81±6.81 


1 


-23.36 


7.11±4.12 


5 


-22.5 


92.4±92.4 


1 


32.7±25.8 


2 


-22.66 


59.9±44.3 


3 


-21.5 
















-20.5 






630±630 


1 


-20.40 


344±344 


1 






Ml4B0 


based on object spectra 








-25.5 


1.83±1.06 


3 


3.05±1.73 


4 


-25.37 


2.49±1.06 


7 


-24.5 


0.63±0.63 


1 


1.19±0.85 


2 


-24.71 


0.93±0.55 


3 


-23.5 


7.66±4.27 


4 


19.6±18.7 


2 


-23.53 


14.2±10.4 


6 


-22.5 


15.0±11.7 


2 


54.6±34.7 


3 


-22.61 


36.6±19.6 


5 


-21.5 


414±414 


1 


630±630 


1 


-21.61 


532±392 


2 



-31 - 



Table 7. Single Power Law Fit 



Dataset 




z,I, i? < 23 


-1.59 ±0.22 
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-2.46 ±0.20 
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Table 8. Double-Power-Law Parameters 
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3000 4000 5000 6000 7000 8000 9000 10000 

Wavelength (A) 

Fig. 1. — Effective tlirougliput curves for the filters used in our quasar color-selection, which include 
the telescope and camera throughputs but do not include atmospheric absorption. The NDWFS 
B\Y^R, and / filter curves are shown in solid blue, orange and red solid lines, respectively. The 
DLS B, V, R (which is identical to the NDWFS R filter), and z filters are shown in blue, green, 
orange and red dashed lines, respectively. We overplot a simulated QSO at z=4.0. 
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Fig. 2. — Positions of the stellar locus color-tracks for the five NDWFS fields (left) and the six DLS 
fields (right). The filled circles indicate the origin of the stellar track and the dashed line traces 
the slope of the locus. Each track is labeled by its corresponding field. 




Fig. 3. — Color-color diagram used in our quasar selection. Small dots are the model quasar colors 
in the NDWFS Bw, R, I filters (left) and the DLS B, R, z' filters (right). Density contours show 
the stellar locus using all stellar objects in each field. The dashed lines show our color criteria for 
selecting quasars in each sample. The open circles show the candidates which passed our candidate 
selection criteria and the filled circles are confirmed quasars. 
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Fig. 4. — Distribution of full width at half-maximum (FWHM) in the i?-band for objects in NDWFS 
J1426-I-3236. We determine the maximum FWHM for selecting stellar sources. Left - FWHM as 
a function of i?-magnitude. The band of stellar sources at FWHM~ 1" is obvious. The dashed 

line shows the maximum FWHM for our candidate selection. Right - Normalized histogram of 
FWHM in two magnitude bins. The dotted line shows objects with R < 22 and the dashed line 
shows objects with 22 < R < 24. We compare these distributions with the FWHM distribution 
of simulated stellar sources (solid-line histogram; see text for a details on the simulations) and 
determine our morphological cutoff at a FWHM that retains 80% of the stellar sources, indicated 
by the vertical solid line at 1.145". 
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Fig. 5. — Contours representing our survey's completeness fraction for the NDWFS survey on the 
left and the DLS survey on the right. Our confirmed quasars arc plotted with filled circles. In most 
of the magnitude and rcdshift range of interest (R < 23, 3.8 < z < 5.2), our completeness is above 
80%, with incompleteness rising at lower redshifts where the Lya emission line is just entering the 
R filter. 
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Fig. 6. — Spectra of DLS and NDWFS quasars, ordered by redshift. The dotted lines show 
expected positions of prominent emission hnes in the ultraviolet: Lya 1216, N V 1240, Si IV 1400, 
C IV 1550, C III] 1909. 
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Fig. 7. — Continued. Spectra of DLS and NDWFS quasars. 
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Fig. 8. — K-correction determined using simulated quasar spectra. The black squares present the 
value of the correction for the NDWFS and DLS i?-band, while the orange circles are the correction 
for the NDWFS I-band and the red triangles are the corrections for the DLS z-band. The ztlo" 
spread is plotted with the dashed curve in the corresponding color. The cyan solid line is the K- 
correction one would obtain using a fixed value of a power-law slope a = —0.5 and no emission line 



or Lyg forest contribution (equation [8]). The blue dashed curve is the K-correction of lRichards et al 



(j2006l ) for the SDSS i'-band. The black diamonds and curve present the K-correction for the SDSS 
i'-band determined using our simulated quasar spectra, with the ±la intervals outlined by the red 
dash-dot curve. As Lya moves through the bandpass, the difference between the correction to the 
R-band and the redder bands increases significantly. 
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Fig. 9. — Absolute magnitude at rest-frame 1450A, M1450, computed using the K-correction derived 
from our library of simulated quasar spectra versus M1450 measured directly from our spectra. DLS 
quasars are marked with asterisks and NDWFS quasars are marked with triangles. The solid line 
is the best fit relationship, with a slope of 1.02 ± 0.09 and an offset of 0.19 it 2.06, suggesting that 
M1450 measured from the spectra are systematically fainter by ~ 0.3 magnitudes. The dotted line 
represents a one-to-one relationship. 
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Fig. 10. — The comoving volume density computed for the quasars in our survey with R < 23. 
Open squares are $(Mi45o, z = 4.15) computed from quasars in the NDWFS survey, open triangles 
are for quasars in the DLS survey, and the filled circles and their corresponding error bars are the 
comoving volume density of the quasars from both surveys combined. The left-hand panel is the 
QLF based on M1450 computed from the z (for DLS) and / (for NDWFS) bands. The right-hand 
panel is the QLF computed using the individual quasar spectra to compute M1450 (as described in 
§5.2p . The dotted line is a single power-law fit, $(M) = ^*10~'^-^(I^~^)^ ^ to the measurements of 
the QLF using the combined survey values. 
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Fig. 11. — The comoving volume density computed for all the quasars in our survey. The symbols 
and lines are the same as in Figure [10] as are the descriptions of the left and right panels. 
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Fig. 12. — The z = 4 .15 quasar luminosity function. Asterisks are the z = 4.25 binned data from 



from 



Fontanot et al. 



Richards et al.l (120061 ) , the open triangles and squares are the combined SDSS and GOODS QLF 



(|2007l ) using at z = 3.75 and z = 4.6, respectively. Our binned QLF data are 
plotted with red circles, where M1450 has been computed from K-corrections to the z and / band 
photometry in the left-hand panel and directly from the spectra in the right-hand panel. The filled 
circles are the mean NDWFS and DLS measurements, while the open circle is the faint bin made 
up of two quasars from the DLS survey. We plot the functional form of the QLF from various 
sources (summarized in the legend). All the points and curves have been evolved to z = 4.15; see 
text for details. 
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Fig. 13. — Ovir selection completeness, fsei, as a function of the quasar's absolute magnitude at 
1450A, Mi45o. The top panel plots the objects with M1450 derived from the K-corrections to z and 
/ magnitudes and the bottom panel plots the objects with M1450 measured from their spectra. The 
horizontal lines are the mean value of the selection function for each luminosity bin. 



